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The braneworld description of our universe entails a large extra dimension and a 
fundamental scale of gravity that might be lower by several orders of magnitude com¬ 
pared to the Planck scale. An interesting consequence of the braneworld scenario is in 
the nature of spherically symmetric vacuum solutions to the brane gravitational field 
equations which could represent black holes with properties quite distinct compared to 
ordinary black holes in 4-dimensions. We discuss certain key features of some braneworld 
black hole geometries. Such black holes are likely to have diverse cosmological and astro- 
physical ramifications. The cosmological evolution of primordial braneworld black holes 
is described highlighting their longevity due to modified evaporation and effective accre¬ 
tion of radiation during the early braneworld high energy era. Observational abundance 
of various evaporation products of the black holes at different eras impose constraints on 
their initial mass fraction. Surviving primordial black holes could be candidates of dark 
matter present in galactic haloes. We discuss gravitational lensing by braneworld black 
holes. Observables related to the relativistic images of strong field gravitational lensing 
could in principle be used to distinguish between different braneworld black hole metrics 
in future observations. 
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1. Introduction 

There has been widespread activity in braneworld gravity in recent times 1 . The 
braneworld scenario of our universe opens up the fascinating possibility of the ex¬ 
istence of large extra spatial dimensions by ensuring that the standard model fields 
are confined to the 3-brane, whereas gravity could also propagate into the higher di¬ 
mensional bulk. Several braneworld models have been studied in the literature, the 
most popular amorng them being the Arkani-Hamed, Dimopoulos and Dvali (ADD) 
model 2 , and the two Randall-Sundrum (RS) models 3,4 . In the ADD model 2 there 
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could be n large compact extra dimensions with radius /, with n > 2 providing a 
possible resolution of the hierarchy problem of particle physics. The RS-I model 5 is 
motivated from similar considerations and consists of two opposite tension branes, 
with our universe stipulated to be the negative tension brane. In the RS-II model 4 
which has been the inspiration for an extensively developed braneworlcl cosmology 5 , 
the AdS5 bulk can have infinite size. Confinement of the standard model fields is 
achieved by a positive tension brane, and a negative cosmological constant for the 
AdS5 bulk with curvature radius l. The resultant modification of the Newtonian 
gravitational potential in braneworlcl models 6 ’' is of the order 1 /r 3 at distances 
r > l where l is the scale of the extra dimension(s). The failure of current exper¬ 
iments using torsion pendulums and mechanical oscillators to observe departures 
from Newtonian gravity at small scales have set the upper limit of l in the sub- 
milimeter region, i.e., l < 0 . 2 m?n 8 . 

It has been recently realized that the possibility of observing signatures of mod¬ 
ified gravity in braneworlcl models with large extra dimensions exists in particle 
accelerator experiments. This is due to the fact that in braneworlcl gravity the 
fundamental 5-dimensional Planck scale M 5 can be much below the 4-dimensional 
Planck scale M 4 . If the scale of the extra dimension is not much below the limit 
obtained from table-top experiments 8 , then the corresponding M 5 could be as high 
as the order of a TeV. Thus it is possible for mini 5-dimensional black holes to 
be produced in particle collisions with centre of mass energy of TeV order 9 . The 
production cross-sections of various higher dimensional black holes in future ac¬ 
celerators such as the LHC has been studied and possible signatures in the form 
of properties of end products of the Hawking evaporation from these black holes 
have been enlisted 10 . High energy cosmic ray showers could similarly produce small 
higher dimensional black holes, the possible signatures from which have also been 
studied 11 . It hence appears that observational signatures of braneworld gravity may 
be a distinct possiblity in the near future from several avenues. 

The analysis of gravitational field equations on the brane is conceptually compli¬ 
cated due to the fact that the propagation of gravity into the bulk does not permit 
the treatment of the brane gravitational field equations as a closed form system 12 . 
This makes the task of studying gravitational collapse on the brane rather difficult 13 . 
Garriga and Tanaka 6 first incorporated the effect of the Kaluza-Klein modes on the 
metric outside a spherically symmetric and static matter distribution on the brane 
in the form of the 1/r 3 correction to the gravitational potential. Since then though 
no exact solution of the full 5-dimensional bulk field equations have been found, 
various solutions representing braneworld black holes have been obtained based on 
different configurations for the projected 5-dimensional Weyl tensor on the brane. 
Among them the one of the Reissner-Nordstrom type with negative tidal charge 14 
originating from the Weyl term has been discussed in details in the literature. The 
projection onto the brane of the 5-dimensional Schwarzschild solution given by the 
Myers-Perry metric 15 could effectively describe a small braneworld black hole of 
size r < l. The properties of such black holes have been investigated in details 16 . 
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The mechanism of Hawking evaporation has been formulated into the brane and 
the bulk as well 17 , and special features have been observed in the interaction of 
5-dimensional black holes with branes 18 . 

Braneworld black holes have been the subject of recent phenomenological in¬ 
terest in the arena of cosmology. Primordial black holes have potentially diverse 
ramifications on several eras of cosmological evolution. If braneworld black holes 
are formed in the early universe, their effect on subsequent dynamics could be not 
only varied but also widely different from those of primordial Schwarzscliilcl black 
holes in standard cosmolgy. This is due to the reasons that braneworld black holes 
have different properties compared to ordinary 4-dimensional black holes, and also 
due to the entirely modified evolution of the very early stages of the universe in 
the braneworld scenario 1,5 . In particular, the Hubble expansion is modified by the 
presence of a term proportional to the square of the energy density on the right 
hand side of the Friedmann equation, which dominates the dynamics during the 
very early high energy phase. Primordial braneworld black holes that could be 
produced due to the collapse of horizon-sized density perturbations, have a lower 
temperature and evaporate slowly compared to the standard 4-dimensional black 
holes as a consequence of their different geometry 19 . Furthermore, the modified 
braneworld cosmological evolution of the universe enables accretion from the sur¬ 
rounding radiation to be effective towards increasing the mass and longevity of the 
black holes 20,21 . Primordial black holes could survive till many different eras in such 
a scenario, thereby contributing to the energy density of the universe. The Hawking 
evaporation products at the end of their life-cycles could have a significant bear¬ 
ing on several cosmolgical processes. Observational results such as the background 
gamma ray spectrum could hence be used to impose constaints on the initial mass 
spectrum of the black holes. These could be again significantly modified 22 com¬ 
pared to those that have been obtained for primordial Schwarzschild black holes in 
standard cosmology. 

The prospect of survival of primordial black holes up to present times in the 
braneworld scenario naturally raises the question as to whether they could be a sig¬ 
nificant fraction of cold dark matter in galactic haloes, and also what role they could 
have in structure formation. In order to address these issues it is first important 
to obtain observational evidences of their existence and to determine their mass 
ranges. The two likely avenues for obtaining observational signatures from black 
holes that may be present in our galactic halo are through gravitational waves from 
coalescing black hole binaries, and through gravitational lensing of light sources 
by the black holes. It has been shown that energy exchange between neighbouring 
black holes that are formed in the high energy braneworld era facilitates the later 
formation of black hole binaries through gravitational interaction 23 . It has been 
also argued that binaries of primordial black holes in the braneworld scenario could 
emit gravitational waves observable by future detectors 24 . The study of gravitational 
lensing by braneworld black holes has been very recently undertaken. The deflec¬ 
tion of light propagating on the brane due to bulk effects has been calculated 25 . 
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The expressions for the weak field bending angle of light in certain braneworlcl met¬ 
rics has been obtained 26 . Further, the various lensing quantities for one possible 
braneworlcl black hole geometry have been obtained and compared to those for the 
Schwarzschild black hole in the weak field limit 27 . The richer phenomenology of 
strong field gravitational lensing such as the positions and magnifications of rel¬ 
ativistic images is being investigated 28 in several braneworld geometries, and the 
values of observational parameters computed for a candidate lense 29 . 

The plan of this review is as follows. In the next section we describe some can¬ 
didate geometries for braneworld black holes. The arbitrariness of the projected 
bulk Weyl term on the brane is responsible for the existence of a number of possi¬ 
ble solutions. We discuss some distinctive properties of a few of them. The aim of 
this article is to highlight the progress made in understanding the impact of hav¬ 
ing braneworld black holes in cosmological and astrophysical processes. In order to 
make this review a bit self-contained we provide a brief description of the essential 
features of the cosmology of the braneworld scenario in section 3. The stage is then 
set for a somewhat detailed analysis of the cosmological evolution of primordial 
braneworld black holes in section 4. Here we try to emphasize the key differences 
from the consequences of a population of primordial Schwarzschild black holes in 
standard cosmology. We begin section 5 with a skeletal description of the theo¬ 
retical framework of gravitational lensing. We then present some recent results on 
braneworld lensing quantities and observables for some of the geometries described 
earlier. The underlying spirit here is the attempt to possibly discriminate between 
different gravity and braneworld gravity models by future observations. A summary 
is presented and some concluding remarks are made in section 6. 


2. Spherically symmetric and static vacuum solutions on the brane 

Obtaining the gravitational field due to a localysed matter distribution on the brane 
has been an involved and challenging task right since the inception of braneworld 
models. This question is the forebearer of the problem of finding the final state 
of gravitational collapse on the brane, which is of central importance concerning 
the existence of black hole solutions in the braneworld scenario. The process of 
gravitational collapse in the braneworld scenario is much complicated compared 
to general relativity because in the former case whereas matter is confined to the 
brane, the gravitational field can also access the extra dimension. The corrections 
to the Newtonian potential of a point mass M at large distances due to the extra 
dimension were calculated to be 3,6 


V(r) 


2 M / 27 2 \ 

M|r \ ^ 3 r 2 ) 


(i) 


The failure of current experiments to detect such corrections at sub-milimeter scales 
have set the upper limit on the curvature radius l of the 5-th dimension as l < 
0 .2m?n 8 . 
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The effect of Kaluza-Klein modes on the metric exterior to a static and spher¬ 
ically symmetric matter distribution on the brane was considered by Garriga and 
Tanaka 6 . They obtained a solution in the weak field limit given by 


dSl 



AMI 2 
3 M|r 3 


dt 2 + 



2 Ml 2 
3M$r 3 


{dr 2 + r 2 d0 2 ) 


( 2 ) 


Note that this solution is quite different from the Schwarzschild metric and that 
the gravitational potential obtained from this metric (2) has 1/r 3 corrections (1) 
compared to the Newtonian potential. Further perturbative studies 7,30 have also 
established that the first weak field correction to the Newtonian potential on the 
brane is proportional to 1/r 3 . 

The projected Weyl term on the brane carries the imprint of Kaluza-Klein 
modes that could be relevant in the process of gravitational collapse. If the Weyl 
term vanishes, then the standard Schwarzschild solution in 4 dimensions can be as¬ 
sumed as the simplest black hole solution on the brane by ‘stacking’ it into the extra 
dimension. Such a vacuum solution of the 4-dinrensional Einstein equation is of the 
‘black string’ type 31 , and can be generalised to the case of a cosmological constant in 
4 dimensions as well 32 . Subsequently, it was shown that the black string is unstable 
to large-scale perturbations 33 . Another solution to the vacuum 4-dimensional Ein¬ 
stein field equations is obtained by setting the 4-dimensional cosmological constant 
to zero (as in Eq.(16), thus obtaining a relation between the brane tension and the 
Ads radius (17)). Since the projected Weyl tensor on the brane is divergence free 
for the vacuum case, one gets for static solutions a closed system of equations given 
by 12 


ft — — E 

K = ° 

= o 


( 3 ) 


Dadhich et al 14 have prescribed the mapping of the 4-dimensional general relativity 
solution with traceless energy momentum tensor of the Einstein-Maxwell type to a 
vacuum braneworld solution in 5 dimensions with the correspondence 


n Cjj u 


( 4 ) 


An exact black hole solution to the effective field equations on the brane of the 
Reissner-Nordstrom type was given with the above correspondence (4) as 14 


dS 2 A 





Q 

r 2 


-1 

dr 2 + r 2 (dfl 2 ) 


( 5 ) 


where Q < 0 is not the electric charge of the conventional Reissner-Nordstrom 
metric, but the negative ‘tidal charge’ arising from the projection on to the brane 
of the gravitational field in the bulk. Since the black hole mass M is the source of 
the bulk Weyl field, the tidal charge Q could be viewed as the reflection back on 
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the brane of the gravitational field of M by the negative AdS5 bulk cosmological 
constant. In the limit r < l, it can be shown that 14,34 


Q = 


Ml 


( 6 ) 


The bulk tidal charge thus strengthens the gravitational field of the black hole. 
It has been further argued 35 that since the back reaction of the bulk onto the 
brane strengthens gravity on the brane, the formation of a black hole as result of 
gravitational collapse is favored as against a naked singularity. The metric with 
negative tidal charge (5) has a spacelike singularity and one horizon given by 


M 





(7) 


which is larger than the Schwarzschild horizon. So the bulk effects are seen to 
increase the entropy and decrease the temperature of the black hole. 

A more general class of spherically symmetric and static solutions to the field 
equations with a 5-dimensional cosmological constant can be derived by considering 
a general line element of the type 


ds 2 = —A{r)di 2 + B(r)dr 2 + r 2 {dn 2 ) 


( 8 ) 


and relaxing the condition A(r ) = i? - 1 (r) used while obtaining the Schwarzschild 
or the Reissner-Nordstrom metrics. Casadio et al 36 obtained two types of solutions 
by fixing either A(r) or B(r), and then demanding the correct 1/r asymptotic 
behaviour for the other in terms of the post Newtonian (PPN) parametrization. In 
the first case, the choice A(r) = 1 — 2 M/(M 2 r) leads to the metric 


ds 4 


-(1 


2 M 

Mir 


)dt 2 


(1 


1 — 3M 
2 Mjr 

2 M \ (i _ M(4/3—1) \ 
Tljr > l^ 4 2 Mfr ) 


+ r 2 (dfl 2 ) 


(9) 


in terms of the PPN parameter (3 which impacts the deflection and time delay of 
light 37 . Note that the above metric was also derived as a possible geometry outside 
a star on the brane 38 . The solution (9) is of the temporal Schwarzschild form having 
a horizon = 2 M/Ml. The corresponding Hawking temperature is given by 36 


= V1-6Q3-1) 
87 tM 


( 10 ) 


Thus, in comparison with Schwarzschild black holes, the black hole (9) will be either 
hotter or colder depending upon the sign of (f3 — 1). 

Alternately, the choice for B(r) of the form R“ 1 (r) = 1 — 2yM/(M|r), in terms 
of the PPN parameter 7 , yeilds the line element 


dS\ = 




dt 2 + 


dr 2 


2 M 
M'lr 


r 2 {dn 2 ) 


7 — 1 + 


2 


(ii) 
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This form of the metric represents a nonsingular wormhole, and has been discussed 
earlier in the literature in the context of 4-dimensional general relativity 39 . Worm- 
hole solutions in the braneworld context have been discussed by Bronnikov et al 40 . 
Furthermore, a class of static, spherically symmetric and non-singular braneworld 
solutions with horizon have been obtained 41 by relaxing the vanishing scalar curva¬ 
ture condition (3) used to obtain the solutions (5),(9) and (11). Stationary solutions 
representing charged rotating black holes have also been found recently 42 . The ar¬ 
bitrariness of the projected bulk Weyl term E and its geometric origin is at the 
root of the variety of braneworld black hole and wormhole solutions since both the 
functions A(r) and B(r) in Eq.(8) have to be determined by it 43 . A specific con¬ 
figuration for the Weyl term with a negative equation of state has been considered 
and the resultant geometry with a singular horizon has been worked out to provide 
one more example of a possible braneworld black hole solution 44 . 

The black hole solution (5) exhibits a 1/r 2 correction to the Newtonian potential 
on the brane in contrast to the weak held correction of 1/r 3 as in the solution (2). 
The solution with tidal charge (5) is reflective of the short distance or strong gravity 
limit where the 1/r 2 correction to the gravitational potential may even dominate 
over the 1/r. This corresponds to the fact that at short distances, braneworld gravity 
is truly 5-dimensional. For short distances r <C l it is natural to consider the 5- 
dimensional Schwarzschild solution as a braneworld black hole candidate given by 15 

dr 2 + r 2 (dn 2 3 ) ( 12 ) 

where the horizon size ?’o is so small (?’o < I) so that the black hole effectively 
“sees” all the spatial dimensions on the same footing. A generalysation to higher 
dimensions 16,45 of the hoop conjecture leads to the above form of the metric as 
a static solution to collapsing matter on the brane. Near the event horizon, the 
black hole would have no way of distinguishing between the bulk dimension and 
the braneworld ones. Numerical simulations for scales sufficiently small compared 
to the AdS scale l seem to also support the existence of static solutions satisfying 
the AdS5 boundary conditions 46 . 

The induced 4-dinrensional metric on the brane near the event horizon of the 
5-dinrensional black hole (12) is obtained by integrating out the extra dimension to 
be 

dS 4 2 = -(l-^)di 2 +(l-^) dr 2 +r 2 (dQ 2 ) (13) 

This 4-dimensional metric is different from the standard 4-dimensional 
Schwarzschild solution as it reflects the 5-dinrensional character of the strong grav¬ 
itational held near the black hole horizon in the form of the 1/r 2 gravitational po¬ 
tential. It is however expected that far from the event horizon the metric (13) would 
approach the standard 4-dimensional Schwarzschild form, as was shown explicitly in 
2 + 1 dimensional braneworld framework 4 '. The properties of small black holes with 
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the geometry given by Eq.(13) have been studied extensively by Argyres et al 16 . In 
general these black holes have lesser temperature and a longer lifetime compared to 
the standard 4-dimensional Schwarzschild black holes. Also since the 5-dimensional 
Planck mass could be much lower compared to the 4-dimensional Planck mass 
(Ms <C M 4 ), these black holes could be produced in particle accelerators 9,10 and 
cosmic ray showers 11 . They could thus provide one avenue of testing higher di¬ 
mensional or braneworlcl physics. Of course, the consequences of a population of 
primordial black holes of the type (13) are potentially rich during many different 
cosmological eras, and these will be described in details in section 4. 

The Myers-Perry black hole 15 in 5 dimensions (12) has been also used in the 
context of braneworlcl models to investigate various effects of its interaction with 
the brane and its radiation onto the bulk. Frolov and Stojkovic 18 have shown that 
a small black hole attached to the brane may leave the brane as the result of a 
recoil due to emission of quanta into the bulk. Such an effect leads to energy loss 
in the brane. This opens up the possibility of observing energy non-conservation 
in particle colliders which may be able to produce these black holes. Radiation by 
rotating 5-dinrensional black holes have also been studied and certain conditions 
have been found when such objects could be stationary 48 . The interaction of 5- 
dimensional black holes with the brane described as a domain wall has interesting 
phenomenological features. In particular, the induced geometry on the brane due 
to a moving bulk black hole has been derived, and an apparent violation of the 
energy condition observed on the brane 25 . Specific features of the interaction of 
rotating black holes with the brane have been studied 49 . Energy flux through the 
horizons of various configurations of the black hole-domain wall system have also 
been investigated 50 . 

Before concluding this section, it needs to be emphasized that although sev¬ 
eral spherically symmetric and static brane black hole solutions with contributions 
from the bulk gravity effects have been found 14,36 ’ 38,41,43 ’ 44 , and further possibili¬ 
ties have been elucidated as belonging to a more general class of black holes 40 , none 
of these are obtained as exact solutions of the full 5-dinrensional bulk field equa¬ 
tions. The analyses of gravitational collapse on the brane have typically yeilded 
non-static solutions 13 . Numerical simulations 46,51 have been inconclusive in this as¬ 
pect. Investigations on stellar metrics 43,52 on the braneworld have been performed 
in order to obtain further insight into the full 5-dinrensional spacetime geometry. 
Understanding the bulk properties have been attempted by extending some partic¬ 
ular braneworld black hole solutions to the bulk 53 . The problem of finding the bulk 
metric which would represent a static and spherically symmetric vacuum solution 
with horizon on the brane remains an open one till date. 


3. Braneworld cosmology 

The cosmology of the RS-II model entails a modified high energy phase in the early 
radiation dominated era of the universe during which the right hand side of the 
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Einstein equation contains a term that is quadratic in the brane energy momen¬ 
tum tensor 1 . Other modifications include the so-called “dark-energy” term which 
is given by the projection of the bulk Weyl tensor. Transition to the standard radi¬ 
ation dominated era takes place when t » t c = 1/2. Such a modified high energy 
evolution has rich consequences for the physics of the early universe 5 . In particular, 
the inflationary scenario is altered, allowing the possibility of steep inflaton poten¬ 
tials to accomplish the desired features. Constraints on the duration of the brane 
dominated high energy phase are enforced by the necessity of conforming to the 
standard cosmological observational features such as nucleosynthesis and density 
perturbations. 

Let us now review briefly some of the essential features of the RS-II braneworld 
cosmology. The effective 4-dimensional Einstein tensor on the brane is given by 1 


G - —t 

M fh 




(14) 


where t m „ is the brane energy-momentum tensor; II M1/ is quadratic in the brane 
EM tensor; and is the projection of the 5-dimensional Weyl tensor. The 4- 
dimensional Planck’s mass AL 4 is related to the gravitational coupling constant n 
and the AdS length l by -jgr = 

For the Friedmann-Robertson-Walker metric on the brane, the Friedmann equa¬ 
tion is given by 


H 2 


87T 

3A/f 




(15) 


with H being the Hubble constant, p the energy density, and k = —1,0,1 represent¬ 
ing open, flat and closed branes, respectively, pkk is the effective energy density 
coming from the bulk Weyl tensor, A = 3M®/4ttM 2 is the brane tension, and A 4 the 
effective 4-dimensional cosmological constant on the brane. The AdS curvature ra¬ 
dius l is given by the bulk cosmological constant A 5 and 5-dinrensional Planck mass 
AL 5 as A 5 = — ( 3 AL|)/( 47 tZ 2 ). The induced 4-dimensional cosmological constant A 4 
is given by 


A 4 




(16) 


Setting A 4 = 0, one obtains a relation between the brane tension and AdS radius 
given by 


A " 1 / 4 = 



(17) 


Nucleosynthesis and CMBR observations constrain the “dark energy” term pkk 
to be negligible compared to the radiation density p 54 . For much of cosmological 
evolution one can neglect pkk, as we will do in the following analysis. 
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Assuming a radiation dominated equation of state, the (k = 0) solutions for the 
Friedmann equation are given by 


= 3 Ml 

32nt(t + t c ) 


(18) 


for the energy density, and 


a = a o 


t(t + t c ) 
to(t o + t c ) 


(19) 


for the scale factor a during the radiation dominated era, and where t c = 1/2 
effectively demarcates the brane dominated “high energy” era from the standard 
radiation dominated era. For times earlier that t c , i.e., t < t c (or p > A), one has 
the non-standard high energy regime during which the radiation density and the 
scale factor evolve as 


and 


3 Mj 
32irt c t 



( 20 ) 


( 21 ) 


respectively. As a consequence, the time-temperature relation also gets modified 
during the brane dominated high energy era, i.e., T oc f^ 1 / 4 . 

On the other hand, if the high energy braneworld regime is matter dominated 
from a time t m , the scale factor grows subsequently like 

/ t \ 1/3 

a = a m I — ) (22) 


But, for times much later than t c , i.e., t » t c (or p « A), one should recover 
back the standard radiation dominated cosmological evolution given by 


PR = 


3 Ml 
327 Tt 2 


and 


a = a o 


{t 0 t c ) 1/2 ) 


1/2 


(23) 


(24) 


The observational success of standard big-bang nucleosynthesis constrains that the 
high energy era be over by the epoch of the synthesis of light elements. However, 
this requirement is satisfied for l < 10 43 / 4 , which is a much weaker bound than that 
obtained from experiments probing the modifications to Newtonian gravity 8 in the 
braneworld scenario. 

Modified expansion in the high energy era has interesting implications for 
inflation 55 . For inflation driven by a scalar field <j> with potential V on the brane, 
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the condition for accelerated expansion of the scale factor (a > 0) is satisfied when 
the equation of state parameter w is 


f 1 + 2p/A \ 
3 \ 1 + p/ A / 


(25) 


In the standard slow role approximation the Hubble rate and the scalar field evolve 
as 


H 2 




V' 
3 H 


(26) 


The braneworld correction term V/2\ enhances the Hubble rate compared to stan¬ 
dard cosmology. This is turn increases friction in the scalar field equation. Thus 
slow roll inflation is favored even for steep potentials 56 . The added advantage of 
such a scenario is that the inflaton field can play the role of quintessence 57 leading 
to a late time acceleration of the universe, as well. Further modifications to the high 
energy expansion in the very early stages can be brought about by the inclusion of 
the Gauss-Bonnet term in the 5-dimensional action 58 . During the very early era the 
Gauss-Bonnet term drives the Hubble rate as H 2 oc p 2 / 3 (Gauss-Bonnett regime), 
which subsequently changes to H 2 oc p 2 (Randall-Sundrum regime), and finally for 
t > t c the H 2 oc p (standard regime) evolution is recovered. Braneworld inflation 
could be accomplished by the Gauss-Bonnet term for very steep potentials such as 
the exponential potential 59 . 


4. Cosmological evolution of black holes 

This section will focus on the consequences of a population of primordial black holes 
on cosmology in the braneworld scenario. The modified features of cosmology during 
the high energy braneworld era have been highlighted in section 2. The black holes 
present in the early universe affect the dynamics of the radiation dominated ex¬ 
pansion through Hawking emission and accretion of the surrounding radiation 20,21 . 
These two competing processes lead to a net energy flow for a single black hole, 
the direction of which determines its longevity. The temperature and the rate of 
Hawking radiation for braneworld black holes are themselves different from those of 
standard 4-dimensional Schwarzschild black holes. The evolution of the individual 
black holes are impacted by the accretion of radiation from the surroundings, which 
is more effective in the braneworld scenario compared to the standard cosmology. 
Since the rate of accretion is governed by the rate of background expansion which 
is much slower in the high energy regime, black holes that are produced earlier 
undergo larger growth. 

The actual rate of accretion and evaporation of course depends upon the partic¬ 
ular geometry of the braneworld black hole. The black holes of interest for cosmo¬ 
logical evolution are produced very early in the universe either due to the collapse 
of overdense regions resulting from inflation generated density perturbations, or 
due to the collision of heavy particles in the primordial plasma. Most of such black 
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holes are expected to be formed with a size small enough (r < l) for the induced 4- 
dimensional Myers-Perry metric (13) to be a good approximation to their geometry 
on the brane. Further, as we are interested in the processes of Hawking evaporation 
and the accretion of radiation, both of which are characterized by the near horizon 
short distance properties of the metric, it may also be pragmatic to consider the 
5-dimensional form of gravity reflected in the near horizon strong field region by the 
geometry in Eq.(13). It is worth noting that the collapse of the ‘tidal charge’ in vac¬ 
uum could also give rise to the same geometry for primordial black holes 14 (M = 0 
in Eq.(5)). For these reasons only this particular form (13) of braneworlcl black hole 
geometry has been considered for the analysis of cosmological evolution 19,20 ’ 21 ’ 22 . 

We first consider the evolution of a single primordial black hole which is formed 
with a sub-horizon mass 19 in the high energy radiation dominated era. Since we are 
considering the 4-dimensional projection of the 5-dimensional Schwarzsclrild metric, 
i.e., 

dSt = - (l- dt 2 + (i_ !^-) 1 dr 2 + r 2 (cm 2 ) (27) 

the horizon radius of such a black hole is proportional to the square root of its mass. 
The mass-radius relationship given by 



which is different from the ordinary 4-dimensional Schwarzschild radius. Various 
properties of such black holes have been elaborated 16 , and the process of Hawking 
evaporation into the bulk and also on the brane have been extensively studied 17 . The 
Hawking evaporation rate, as for the case of standard black holes, is proportional to 
the surface area times the fourth power of temperature. The Hawking temperature 
is given by 


Tbh 


1 

27 T tbh 


(29) 


Therefore, as a consequence of the mass-radius relationship (28), such black holes 
are colder and long-lived compared to 4-dimensional Schwarzschild black holes. 

A black hole formed in the early radiation dominated era of the universe accretes 
the surrounding radiation. In standard cosmology the effectiveness of accretion in 
the growth of black hole mass is restricted because of the fact that the mass of the 
individual black holes could grow at nearly the same rate as that of the cosmolog¬ 
ical Hubble mass Mh. i.e., Mh ~ M ~ t. Black holes produced in the radiation 
dominated era cannot be formed with a size much smaller than the Hubble radius, 
since otherwise pressure forces could hinder the collapse process. However, in the 
radiation dominated high energy phase of the braneworld scenario, the Hubble mass 
grows as Mh ~ t 2 , whereas the growth of black hole mass due to accretion is given 
by M ~ t B (with B < 2 /pi) 20,21 . Hence, sufficient energy is available within the 
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Hubble volume for a black hole to accrete in the high energy braneworld regime. 
The exact efficiency of accretion though depends upon complex physical processes 
involving the mean free paths of the particles comprising the radiation background 
and the thermal properties of the radiation in the non-trivial geometry near the 
event horizon 21 . Any peculiar velocity of the black hole with respect to the cos¬ 
mic frame further impacts the rate of accretion. In the absence of a universally 
accepted approach of determining the precise accretion rate, it is usually taken to 
be as equal to the product of the surface area of the black hole, the energy density 
of the background radiation, and an efficiency factor ranging between 0 and l 21 . 

Taking into account these effects of accretion and evaporation together, the rate 
of change of mass M of a braneworld black hole is given by 

M = 4tt r 2 BH (-gbraneCrTgn + fpn) (30) 


where gbrane is effective number of particles that can be emitted by the black hole 
(we assume that the black holes can emit massless particles only and take gbrane = 
7.25 19 ), f is the accretion efficiency (0 < / < l) 21 , and a is the Stefan-Boltzmann 
constant. The black hole also evaporates into the bulk, with a rate proportional 
to 4:TTrg H gbuikTg H . However, this term is subdominant even for very small black 
holes 19 , and has negligible effect on their lifetimes. 

Substituting the expressions for the black hole radius (Eq.(28)), the temperature- 
radius relation (29), and the energy density of radiation (20), the black hole rate 
equation (30) in the radiation dominated high energy braneworld era can be written 


M = - AM2 ‘ 


BM 


Mt c t 

where A and B are dimensionless numbers given by 

3 


A ~ 
B ~ 


(16) 3 7t 

2 / 

7r 


The exact solution for the black hole rate equation is given by 


20 


M(t) = 


M 2 - 


2 AMl t 0 
2 B - 1 T c 



S 2 B 1/2 

t \ 2 AM? t 


2B - 1 t c 


(31) 


(32) 

(33) 


(34) 


with M 0 being the formation mass of the black hole at time t = t 0 . If the black hole 
is formed out of the collapse of horizon or sub-horizon sized density perturbations, 
the formation time and mass are related by 19 


— 

U 


1 Mo 
4 1 M 4 


1/2 


l 

k 


1/2 


(35) 
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It has been argued 20 that a black hole so formed continues to grow in size by the 
accretion of radiation during the high energy radiation dominated era, with its mass 
increasing as 


M(t) /A 

Mo \to) 


(36) 


This result is however sensitive to the accretion efficiency, since a more careful 
analysis 21 shows that for / < 7r/4, the growth due to accretion could come to a halt 
during the high energy regime itself. 

Within the context of standard cosmology, the evolution of a population of pri¬ 
mordial black holes exchanging energy with the surrounding radiation by accretion 
and evaporation has been studied by several authors 60,61 . The analysis of this prob¬ 
lem is simplified by assuming that all the black holes are formed with an average 
initial mass M$ at a time to when the fraction of the total energy density in black 
holes is /3bh , and the number density of black holes is With these assump¬ 

tions, the coupled cosmological equations for the radiation density pn(t), the matter 
density in the black holes and the scale factor a(t ) are integrated to 

give the complete cosmological evolution. Note however, that the black holes are 
produced with an initial mass spectrum in any realistic scenario of black hole for¬ 
mation in the standard cosmology 62 , and it is expected that the same would be the 
case in the braneworld scenario as well. The effect of a mass distribution can be 
incorporated into the cosmological evolution in the standard scenario by introduc¬ 
ing an additional differential equation for the distribution function and specifying 
further initial conditions related to it 61 . However, since not much is known about 
the formation processes of black holes in braneworld cosmology, the study of only 
a few basic features of their cosmological evolution under simplifying assumptions 
has been undertaken in the literature till date. 

The number density of black holes nsHit) scales as a(t )~ 3 , and thus for a radi¬ 
ation dominated evolution on the brane, one gets 

( nBH(t)/riBH{to )) = (t 0 /t) 3/4 (37) 


since a(t) <x f 1 / 4 . The net energy in black holes grows since accretion dominates 
over evaporation. The condition for the universe to remain radiation dominated 
(i.e., pBH(t) < pR{t )) at any instant t can be derived to be 20 


(3bh < 


(fo/i ) B+1/4 
1 + (foA) B+1/4 


(38) 


If the value of /3bh exceeds the above bound, there ensues an era of matter (black 
hole) domination in the high energy braneworld phase. Such a phase of matter dom¬ 
ination should definitely be over by the time of nucleosynthesis for the cosmology 
to be viable. 
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Let us first describe the situation when the cosmology stays radiation dominated 
up to the time when brane effects are important, i.e., t < t c . From Eq.(38), this 
requires 


/ \ B+l/4 

Pbh / to | 

1 — PbH Y tc J 


(39) 


Further, the black holes should remain small enough, i.e., (M/M 4 ) < ( 37 r/ 4 )(i/i 4 ) 
for the 5-dimensional evaporation law to be valid 19 . These criteria can be used to 
put an upper bound to the average initial mass 20 


Mo ( 3tt \ 2 ~ B t c 
Mi y 4(2 v / 2) b J U 


(40) 


The growth of the black holes in the radiation dominated era due to accretion 
slows down with time, since the surrounding radiation density gets diluted. The rate 
of evaporation is also insignificant for a wide range of at this stage since the 
black hole masses could have grown by several orders of magnitude from their initial 
values. There ensues an era during which the black hole mass stays nearly constant 
over a period of time, as is the case for standard cosmology 60 . The accretion rate is 
smaller for the braneworld case since the surface area is proportional to M instead 
of M 2 for 4-dimensional black holes. Moreover, the evaporation rate is (oc Af _1 ) 
instead of M~ 2 . Hence, the black hole mass will stay for while near a maximum 
value M max reached at time before evaporation starts dominating. The expression 
for the lifetime t en d of a black hole in this scenario is given by 20 


tend ^ f 

ti A l M 4 J ti l t c ti I 


(41) 


It is important to note that the modified evaporation law also contributes to the 
increased lifetime for braneworld black holes 19 . However, the effect of accretion 
is more significant, as can be seen by comparing the lifetime of a 5-dimensional 
black hole in the presence of accretion to the the lifetime due to purely an altered 
geometry 20 . Depending upon the values taken by the parameters t c and Z, one could 
obtain several interesting examples of primordial black holes surviving till various 
cosmologically interesting eras 20,21 . One particular choice worth mentioning is that 
of a black hole formed with an initial mass M 0 = 10 8 M 4 ~ 10 3 g, that will survive 
up to the present era if one chooses (l/U) — 10 30 . 

Further interesting phenomena occur in the interaction of two or more neigh¬ 
bouring black holes mediated by the surrounding radiation in the radiation dom¬ 
inated high energy era. These black holes exchange energy via the processes of 
evaporation and accretion with the radiation bath, and through it, with each other. 
The evolution equation for a black hole gets modified due to the presence of one or 
more neighbouring black hole(s) with the addition of an extra source (sink) term due 
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to the evaporating (accreting) neighbour(s) over the average radiation background. 
The black hole equation (31) now becomes 

Bm.j 1 

mi 


mi 


' 


~~9- 


fi/ 2 


(42) 


where t = Am\t/t c , and 


9 = 


4A 1 / 2 / 1 4 
3ir \do 


to\ 

Uj 


1/2 


( — \ 

\ ^4 / 


1/2 


(43) 


can be dubbed as the ‘coupling’ parameter between two black holes. The physical 
distance between two such neighbours increases initially with the Hubble expansion. 
Forming out of horizon collapse, the initial mass ratio of two such black holes is 
on average proportional to the square of the ratio of their formation times (from 
Eq.(35)). By sudying the effect of the interaction term it is possible to show that 
an initial mass difference between two such black holes can never decrease during 
the radiation dominated era 23 . 

The energy exchange between the black holes and the surrounding radiation 
always causes mass disequilibration between neighbours. Such mass differences fa¬ 
cilitate the formation of binaries during the standard low energy phase via three- 
body gravitational interactions. A formation mechanism for primordial black hole 
binaries in standard cosmology has been investigated 63 which shows that if equal 
mass primordial black holes are inhomogeneously distributed in space, three-body 
gravitational interactions could lead to the formation of binaries. The scheme of 
binary formation in the braneworlcl scenario 23 could become more effective when 
combined with the scheme based on spatial inhomogeneities discussed in the con¬ 
text of standard cosmology 63 . Coalescing braneworld black hole binaries may emit 
gravitational waves amenable for detection by the next generation detectors 24 . 

Let us now describe the case of matter (black hole) domination in the high 
energy phase. The onset of such a matter dominated era is derived to be 20 


theq 

to 


1 — Pbh 
Pbh 


= 7 B 


(44) 


The mass of a black hole at th eq is given by M(th eq /M 0 ) = For t > t heq the 
Hubble expansion is essentially driven by the black holes (p = 0) which dominate 
over radiation. Since the number density of black holes scales as matter (n bh (t) oc 
a -3 ), for t <t c one has H oc pbh , and thus the scale factor grows as 

a(t) ~ t 1/3 (45) 


During this era, the radiation density pn is governed by the equation 

^Pi{(f)a 4 (f)^ = (46) 

where the contribution from accreting black holes is comparable to the normal 
redshiting term {p _r ~ a~ 4 ) because at this stage the black holes dominate the 
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total energy density. Some further analysis leads to the following expression for the 
radiation density 20 


PR{t) ~ 1b 1 o) 



PbhB 
B+ 1/3 


7s /4 p(*o) 



1 — B 


(47) 


The black hole mass grows as 


M(t) = Moyfexp 


C 

3B + - 


7 r /4 


- 3 B 



1/3 


1/4 
7B 


C 



(48) 


where C = . As in the case of standard cosmology, the accretion regime lasts 

for a brief duration in the matter dominated phase, beyond which the black hole 
evaporation starts to play a significant role. 

The universe gets reheated as pnit) increases with time. The stage of black hole 
domination lasts up to a time t r ( pB.{t r ) = riBH{t r )M(t r )). Subsequently, radiation 
domination takes over once again. One can derive 20 

h. ~_ - _I_ — (M±\ (49) 

tr 25 b 7 b \yM 0 J 


with 


. _/B + l/3\ %H+ 

■ 5fl = (TT3^J 


(50) 


A stringent restriction is imposed by demanding that the standard low energy cos¬ 
mology (for t > t c ) should emerge as radiation dominated. By requiring that the 
era of black hole domination be over before t c , he., t r < t c , one gets a lower bound 
on /3bh from Eq.(35), i.e., 


(3bh > 


!; <b+1/3) 


3 

3B + 1 


B+ 1/4 


(51) 


The evaporation time of black holes in this scenario has also been calculated 20 . The 
black hole lifetime is given by 


tend ^ I A/o \ t c 2B 

~u M4 u 1b 


(52) 


The effect of accretion is less significant in prolonging black hole lifetimes in this 
case, as borne out by the following example. For instance, taking the value of the 
AdS radius to be l/l 4 ~ 10 20 , and (3bh = 10“ 3 , black holes with Mo ~ 10 12 g 
evaporate during the present era. This is only to be expected, since early mat¬ 
ter domination results in larger Hubble expansion rate which further restricts the 
availability of radiation for the black holes to accrete. 
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It thus turns out that it is possible to have primordial black holes formed in 
the high energy era of the braneworld scenario to survive up to several cosmologi- 
cally interesting eras. Accretion of radiation during the radiation dominated era is 
primarily responsible for the increased longevity of braneworld black holes, though 
their 5-dimensional geometry also contributes to a slower rate of Hawking evapo¬ 
ration. It is worth noting here that primordial black holes in certain models could 
also accrete the energy of a cosmological scalar field 64 . Such an effect could lead 
to further growth of these black holes beyond the early radiation dominated era, 
thus pushing up their lifetimes further. The implications for a population of pri¬ 
mordial black holes in cosmology are diversely manifold. At any particular era, the 
surviving black holes would contribute a portion to the the total energy density as 
dark matter in the universe. If the black holes are produced with an initial mass 
spectrum, then one would have evaporating black holes at different eras. Hawking 
radiation from these evaporating black holes would on one hand produce all kinds 
of particles including heavier ones which could lead to baryogenesis 60 , and on the 
other, could contribute significantly to the background photons, thus diluting the 
baryon to photon ratio. 

Observational constraints impacting different cosmological eras could be used 
to impose restrictions on the initial mass spectrum of braneworld black holes in 
a manner similar to the primordial black holes in standard cosmology. Clancy et 
al 22 have shown how standard constraints are modified in the case of braneworld 
cosmology. To simplify the treatment, one can assume that radiation domination 
persists up to t c , and that the accretion of radiation is possible only up to t c . A black 
hole mass fraction ajf in terms of the radiation density can be defined (related to 
the mass fraction /3bh) as 

u\ PBHit) 

= - 

Pt 

“Mo (to) = 1 3 p BH (53) 

Similarly, a ‘final’ mass fraction a(t e vap) is defined at the end of the black hole 
life-cycle, since the black holes radiate most of their energy towards the very end 
of their lifetimes. Observational constraints on a(f e vap) are considered at different 
cosmological epochs, given by 


a(t) < L iD (t) (54) 

for standard 4-dimensional black holes, or 

a(t) < L 5D (t) (55) 


for braneworld black holes. Then these constraints are evolved backwards to con¬ 
strain initial mass spectrum. Accretion in the high-energy phase leads to a(t) oc 
MBH{t)a(t). Therefore the constraints on the initial black hole mass fraction are 
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given by 

L 4 d = a i < [ ]4_D-b4£>(f) (56) 

a(t.) 

for standard cosmological evolution, and 

L 5D =on< [-^}5DL 5D (t) (57) 

in the braneworld scenario. Any astrophysical or cosmological process to be con¬ 
strained at certain epoch is dominantly affected by the PBHs with lifetimes of that 
epoch. The constraints on primordial black holes in standard cosmology thus get 
modified to 22 

5-16 B 

^5D __ A5D(tevap) / l \ 16 ~ 8Jj (gg) 

^4 D -^4D(^evap) \^min/ 

1 /3 

where Z m ; n oc t e vap- The departure from standard constraints is sensitive to the 
accretion efficiency, as is expected. 

A more detailed study on how the initial mass spectrum of the black holes is 
distorted due to braneworld accretion has been undertaken by Sendouda et al 65 . The 
diffuse photon background emitted by the spectrum of black holes has been shown 
to be modified in accordance with the mass spectrum of the black holes. These 
results have been compared to the observed diffuse photon background to obtain 
bounds on the initial black hole mass fraction, the scale of the extra dimension, and 
the accretion efficiency by these authors. The observed number density of massive 
particles could also be used to obtain bounds on the initial mass fraction of the black 
holes, as in the case of Scliwarzschild primordial black holes in standard cosmology. 
Further constraints on the scale of the extra dimension have been derived 66 by 
considering the recent observation of sub-Gev galactic antiprotons 67 as originating 
from braneworld black holes present in our galaxy. A relevant issue is to investigate 
whether primordial braneworld black holes could contribute to a significant fraction 
of cold dark matter. If so, the direct searches of cold dark matter compact objects 
through gravitational lensing might be able to reveal their presence in galactic 
haloes. In the next section we will report on the specific features of the analysis of 
gravitational lensing for several braneworld black hole metrics. 

5. Gravitational lensing by braneworld black holes 

The braneworld scenario implies the modification of Einstein’s general relativity at 
short distances or strong gravitational fields. The bending of light due to the grav¬ 
itational potential of a massive object is one of the first predictions of the general 
theory of relativity. Its application in the phenomenon of gravitational lensing 68 has 
potentially diverse possibilities. Gravitational lensing in the weak field limit 69 is till 
date one of the most widely used tool in observational astrophysics and cosmology. 
On the other hand, strong field gravitational lensing, though limited in observa¬ 
tional utility because of presently inadequate instruments, remains our ultimate 
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scope for exploring the physics of strong gravitational fields. The general technique 
for analysing strong gravitational lensing for spherically symmetric metrics has been 
developed by Bozza 70,71 who has also formulated useful connections between ob¬ 
servational quantities like fluxes and resolutions and the metric parameters. Strong 
gravitational lensing is also endowed with richer phenomenological features like 
relativistic images 72,73 and retrolensing 74,75,76 . It would be thus worthwhile to in¬ 
vestigate if the results of strong gravitational lensing could be used for probing the 
modifications to general relativity made in braneworlcl geometries. Such studies are 
also motivated from the possibility of producing braneworld black holes in future 
accelerators 9,10 . 

Analysis of the trajectories of light and massive particles in various braneworld 
and higher dimensional metrics have been undertaken recently. Kar and Sinha 26 
obtained the bending angle of light for several brane and bulk geometries. Frolov et 
al 78 found certain non-trivial features about the propagation of light in the Myers- 
Perry 15 metric for a 5-dimensional black hole solution. This solution represents 
primordial black holes that could be produced with size r < l in the early high 
energy era of the RS-II braneworld scenario 19 . It was shown that such black holes 
could grow in size due to accretion of radiation, and consequently survive till much 
later stages in the evolution of the universe 20 . With a suitable choice of parame¬ 
ters, some of these black holes could also exist in the form of coalescing binaries in 
galactic haloes at present times 23 . The weak field limit of gravitational lensing was 
studied for the Myers-Perry metric and certain notable differences from the stan¬ 
dard Schwarzschild lensing were found 2 '. Thereafter, Eiroa 28 analysed the strong 
field lensing and retrolensing effects for the Myers-Perry black hole. Strong field 
gravitational lensing in a couple of other braneworld metrics has been discussed by 
Whisker 29 and some lensing observables have been computed using parameters for 
the galactic centre black hole. 

For a spherically symmetric metric 

ds 2 = — A(r)dt 2 + B(r)dx 2 + C(r) (dfi 2 ) (59) 


where the asymptotic forms of the functions A(r) and B{r) have the standard 1/r 
form, and C(r ) —> r 2 asymptotically, the general formalism of strong field gravita¬ 
tional lensing has been worked out by Bozza 71 . It is required that the equation 

C(r) A{r) 1 ’ 

admits at least one positive solution the largest of which is defined to be the photon 
sphere r ph . 

A photon emanating from a distant source and having an impact parameter u 
will approach near the black hole at a minimum distance ?’o before emerging in a 
different direction (see Figure 1). The closest approach distance is given in terms of 
the impact parameter by 



u = 


(61) 
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where the functions C and A are evaluated at ro. The deflection angle of the photon 
in terms of the distance of closest approach is 


a{ro) = I(r 0 ) - 7r 


,(r " ) = / 

^0 


2 yfBdr 



_C_Aq. _ 1 

C 0 A 1 


(62) 

(63) 


The weak field limit is obtained by expanding the integrand in Eq.(63) to the first 
order in the gravitational potential. This limit however is not a good approximation 
when there is a significant difference between the impact parameter u and the 
distance of closest approach ro, which occurs when A(tq) significantly differs from 
1, or C(r o) from r p . By decreasing the impact parameter, and consequently the 
distance of closest approach, the deflection angle increases beyond 27 t at some stage 
resulting in one or more photonic loops around the black hole before emergence. 
Further decrease of the impact parameter to a minimum value u m corresponding to 
the distance of closest approach r 0 = r p h results in the divergence of the deflection 
angle (integral in Eq.(63)), which means that the photon is captured by the black 
hole. Strong field gravitational lensing is useful for studying the deflection of light 
in a region starting from just beyond the photon sphere up to the distance where 
the weak field approximation approaches validity. 

In the general analysis of strong field gravitational lensing it becomes necessary 
to extract out the divergent part of the deflection angle. In order to do so two new 
variables y and z are defined as 71 


y = A(r) 

y-yo 

z = 1 - 

1 - 2/0 


(64) 

(65) 


where yo = Aq. In terms of these variables, the integral (63) in the deflection angle 
is given by 


1 



R(z,r 0 )f{z,r 0 )dz 


R(z,r 0 ) 

f(z,r 0 ) 


2 VBy 

CA! 


(1 - 2/o) 


1 


\Jyo- [(l - 2/0) 2 + yo] %- 


( 66 ) 

(67) 

( 68 ) 


where all functions without the subscript 0 are evaluated at r = 
A _1 [(1 — yo) z + yo]- The function R(z,ro) is regular for all values of z and ro, 
while /(z, r 0 ) diverges for z —> 0. 

The deflection angle can be written as a function of 9 = u/Dd , where 6 is the 
angular separation of the image from the lens, and Dd is the distance between the 
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lense and the observer (See Figure 1). In order to do so, an integral 


1 

b R = J g(z,r m )dz 


(69) 


is defined, where 

g(z, Tm) = R{z, r m )f(z, r m ) - R(0, r m )f 0 (z, r m ) (70) 

The expression for the deflection angle is given as a function of 9 = u/Dd as' 1 

'0D d 


a(6) = —a log 

_ _ R( 0, r m ) 
a ~ 2 ^ 

b = — 7r + b R + a log 


- 1 


2ft, 

Um 


(71) 

(72) 

(73) 


where R(0,r m ) and b R are given be Eqs.(67) and (69), respectively, and /3 m is 
defined as 


Pm = 


C m (1 - Hmf (C'^y m - C m A" (r m )) 


2 VlC'm 


(74) 


In strong lensing there may exist n relativistic images given by the number of 
times a light ray loops around the black hole. The positions of these are obtained 
as solutions of the lense equation given by 


t anS = t and - -^-[tanO + t an(a — 0)] 


(75) 


for specific positions of the source and the lense respective to the observer, and using 
the the value of the deflection angle from Eq.(71). The relativistic images formed 
by light rays winding around the black hole are highly demagnified compared to the 


S I 



Fig. 1. Gravitationl lensing for point like mass object M. A light ray from the source S passes 
the lens with an impact parameter u , and is deflected by an angle a. The observer sees an image 
I of the source at the angular position 6. 
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weak field images. When the source, the lense and the observer are highly aligned, 
it is possible to obtain the most prominent of the relativistic images' 2 . Hence, the 
analysis of strong lensing is usually restricted to the case when both 5 and 6 are 
small 71 , though the general case for arbitrary positions can also be analysed 76 . With 
the above restriction on the values of 5 and 9 , a light ray will reach the observer 
after winding around the lense n number of times only if the deflection angle a is 
very close to a multiple of 27 t. Substituting a = 2nn + A a n in Eq.(75), one gets 

S = 9- ^Aa„ (76) 

1C S 

The position of the n-th relativistic image 9 n can hence be obtained as a solution 
of the lense equation (76) as 71 


ttm^n ( $ 


e n = ^(l + e n 

J->d 


,(l+e 


D d 


D„ 


uD ds D d 

where u m is the minimum impact parameter, and e n is given by 

p — Jb-2nir)/a 
c-n — ' 

The magniheation ji n of the n-th relativistic image is given by' 1 


Idn — 


(5/9)8889 


\e n - 


2 ^n( 1 T &n)D s 
a5D ds D 2 , 


(77) 


(78) 


(79) 


The above formula for magniheation is valid under the approximation of a point 
source. However, for an extended source the magniheation at the image positon can 
also be derived by integrating over the luminosity profile of the source 75 . 

It is useful to obtain the expressions for the various lensing observables in terms 
of the metric parameters. For n —> oo an observable 9^ can be defined 71 represent¬ 
ing the asymptotic position approached by a set of images. The minimum impact 
parameter can then be obtained as 


u m = D d 9oo (80) 

In the simplest situation where only the outermost image 9\ is resolved as a sin¬ 
gle image, while all the remaining ones are packed together at 9 oo, two lensing 
observables can be dehned as 71 


5 = - 0oo (81) 

representing the separation between the hrst image and the others, and 

n = -5^- (82) 

E 9-n 
n =2 

corresponding to the ratio between the hux of the hrst image and the flux coming 
from all the other images. 



February 2, 2008 17:6 WSPC/INSTRUCTION FILE review 


24 A. S. Majumdar, N. Mukherjee 


In terms of the deflection angle parameters a and b , these observables can be 
written as 71 


5 = 0 oo e b / a_27r / a 

(83) 

(M 

0 ) 

II 

(84) 


The above equations (83) and (84) can be inverted to express a and b in terms of the 
image separation S and the flux ratio 1Z. Therefore the knowledge of these two ob¬ 
servables can be used to reconstruct the deflection angle given by Eq.(71). The aim 
of strong field gravitational lensing is to detect the relativistic images corresponding 
to specific lensing candidates and measure their separations and flux ratios. Once 
this is accomplished, the observed data could be compared with the theoretical 
coefficients obtained using various metrics. A precise set of observational data for 
strong gravitational lensing, if obtained, could therefore be able discriminate be¬ 
tween different models of gravity. In the braneworlcl scenario the computation of 
the above parameters has been performed for two metrics 29 taking the black hole 
at the centre of our galaxy as a potential candidate. 

We will now consider examples of the various lensing quantities defined above 
for some of the possible braneworlcl black hole geometries discussed in section 2 . 
It is instructive to compare the braneworld lensing quantities with the standard 
Schwarzschild ones which for strong gravitational lensing are given as follows 71 . 
Choosing the Schwarzschild radius r s = 2 M /Mf as the unit of distance, the photon 
sphere is given by 

r P h = | (85) 


for Schwarzschild lensing. The corresponding minimum impact parameter is 


— 


3-\/3 


( 86 ) 


The coefficients a and b defined in Eqs.(72) and (73) are given by 


a = 1 

b = — 7 t + 2log[6(2 — \/3)] + log[ 6 ] 


(87) 

( 88 ) 


The deflection angle (71) is obtained in terms of the parameters a and b to be ' 1 


a(0) = -log 


(29D d 

UV3 


■ lop[216(7 - 4-v/3)] - 7 T 


(89) 


In the weak field limit, the expression for the bending angle is given by 

AM 

Q^vreak 


Mb o 


(90) 
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For the analysis of gravitational lensing by braneworld metrics, let us first con- 
sider the Garriga-Tanaka weak field solution 6 given by Eq.(2) in isotropic coordi¬ 
nates. The metric in terms of the standard coordinates can be written as 


ds\ 



AMI 2 \ 

3Mjr 3 ) 


dt 2 + 



1 + 


2M 

mJ7 


Ml 2 
3 M-fr 3 
2MP 
3 MjV 3 


-2 

dr 2 + r 2 (dfl 2 ) 


(91) 


The formal expression for the radius of the photon sphere r p / 1 can be obtained as a 
function of M and l using Eq.(60) to be 

r s 3 1 / 3 r 2 

7ph = + 2(3r 3 - 20 l 2 r s + 2y/Wy/Wl 4 T* - 3Z 2 r 4 ) 4 / 3 

, (3r s 3 - 20 l 2 r s + 2 v / 10 \/ 10 / 4 rf - SPri) 1 ' 3 
+ - 2 ^- (92) 

where r s = 2 M/M 2 . However, it can be seen from Eq.(92) that no real solution for 
r p h exists for admissible values of l and r s . This is to be expected since the metric 
(2) represents a weak field solution. The weak field limit of the bending angle was 
obtained to be 26 


AM AMI 2 

owfc - ^ (93) 

where ro in this case is the isotropic coordinate equivalent of the distance of closest 
approach ro in standard coordinates. 

The analysis of strong field gravitational lensing can however be performed in 
other braneworld geometries. For example, let us consider lensing in the metric with 
tidal charge (5) given by 14 


dS\ 






dr 2 + r 2 (dfl 2 ) 


(94) 


which resembles the Reissner-Nordstrom metric, but with Q < 0 in the braneworld 
context. Again using units of distance 2M/M 2 , it is straightforward to obtain the 
expressions for the photon sphere and the minimum impact parameter given by 71 


(3 + \/9 - 32Q) 

4 

(3 + V9-32Q) 2 
4^x73 -8 Q + 79 ~ 32Q 


(95) 


(96) 
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The coefficients a and b in the deflection angle are given by 71 


h \/r r ,h 2Q 


b = —7 r + 21o</[6(2 — 73)] + 


^/(3 r ph) r ph 9Q r ph + 
8Q (73 — 4 + \og[6(2 — 73 


(r ph - Q) 2 [(3 - T p h)r ph - 9 Qr ph + 8<5 2 ]alo^[2] 
( r P h ~ 2Q) 3 (r 2 /l - r ph + Q) 


( 97 ) 


(98) 


In terms of the above coefficients one obtains the complete expression for the de¬ 
flection angle using Eq.(71). The weak field limit of the bending angle was derived 
to be 26 


^weafc — 


1 

m! 


37 tQ \ 


4 M 


16M?’o J ro 

Note that the bending angle is always positive because of negative tidal charge Q 
unlike the electric charge of the Reissner-Noredstrom metric. 

We next consider the braneworld solution (9) in terms of the PPN parameter /3 
given by 36 

1 _ 3 M 

)dt 2 + - — 1>N + r 2 (dfl 2 ) (100) 


ds\ = -(1 - 


2 M 

Mir 




Af (4/3-1) 
2 Mh- 


For the above metric the expressions for the radius of the photon sphere r p h and the 
minimum impact parameter are of course similar to those for Schwarzschild lensing 
given by Eqs.(85) and (86), as is easy to see using Eqs.(59), (60) and (61). The 
expression for the deflection angle can be derived in terms of the coefficients a and 
b which for the metric (100) are given by (setting the unit of distance as 2M/M 2 ) 


V3 

\/6-(4/3-1) 


( 101 ) 


b = — 7T - 


2V3 


\/6 - (4/3 - 1) 


%[6(2 - V3)] 


73 


76 - ( 4/3 - 1 ) 


lop [6] 


( 102 ) 


Kar and Sinha 26 obtained the weak field limit of the bending angle for the metric 
(100) to be 


^wea/c 


2M(l + /3) 
ro 


(103) 


Note that the standard Schwarzschild expressions are recovered in Eqs.(lOl), (102) 
and (103), as should be, for (3 = 1. 

The lensing observables S and 1Z corresponding to the ratio between the flux of 
the first image and the flux coming from all the other images as defined in Eqs.(83) 
and (84) respectively, can be computed using the expressions for u m , a, and b for 
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particular geometries. One can obtain the magnitudes of S and 7Z for particular 
lensing candidates using the known values for their masses and distances. For the 
black hole located at the centre of our galaxy at a distance of Dd = 8.5kpc and with 
mass M = 2.8 x 10 6 M 0 , the position of relativistic images in Schwarzschild strong 
lensing was first computed by Virbhadra and Ellis 72 , and the observable parameters 
S and 7Z were computed by Bozza 71 . Whisker 29 has computed using the above black 
hole as a candidate lense the values for these observables for two possible braneworld 
gemetries given by the metric (94) with tidal charge 14 and another solution 44 . 


Table 1. Estimates for the strong field lensing angle coefficients and observables for the black hole at the 
center of our galaxy with mass M = 2.8 x 10 6 Mq and Dd = 8.5 kpc, for standard Schwarzschild and two 
different braneworld geometries. (r m = 2.5\og7Z.) 


Observables 

Schwarzschild 

metric 

Brane metric 
with tidal charge Q 

Brane metric 

with PPN parameter (3 

Q = -0.1 

Q = -0.2 

Q = -0.3 

p = 1 + lO -4 

1 

o 

1 

II 

#oo(m arc sec) 

16.87 

17.87 

18.92 

19.65 

16.87 

16.87 

S(/j, arc sec) 

0.0211 

0.0142 

0.0102 

0.097 

0.02115 

0.01923 

r m (magnitudes) 

6.82 

7.02 

7.2 

7.37 

6.818 

6.887 

Urn ft's 

2.6 

2.75 

2.9 

3 

2.6 

2.6 

a 

1 

0.9708 

0.938 

0.925 

1.00006 

0.9999 

b 

-0.4002 

-0.612 

-0.747 

-0.819 

-0.402 

-0.429 


In Table 1 we display the values of these lensing quantities 9^, <S, TZ, u m /r s , a, 
and b first for Schwarzschild lensing. We compare the values of the lensing ob¬ 
servables with those obtained for lensing with two braneworld geometries, i.e., the 
metric (94) with tidal charge 14 , and the metric (100) with the PPN parameter /3 36 . 
We choose three negative values of the tidal charge Q. The two values of the pa¬ 
rameter (3 on either side of the Schwarzschild value (/3 = 1) are chosen to maintain 
observational compatibility with the solar system tests of the Nortvedt effect 3 '. The 
measurement of the observable 9 oo involves a microsecond resolution which should 
in principle be attainable by the very long baseline interferometry projects such as 
MAXIM 77 . Though the actual identification of faint relativistic images would be 
extremely difficult in practice due to the inherent disturbances 72 , an accurate mea¬ 
surement of 9oo would be able to distinguish the Schwarzschild geometry from the 
braneworld RN-type one. However, to unambiguously determine the exact nature 
of the black hole through the lensing angle coefficients o and 6, one has to measure 
the observables S and TZ. Since this involves the resolution of two faint images sep¬ 
arated by ~ 0.02/z arc sec, such an observation would need a leap of technological 
development over the present astronomical facilities 29 ’' 1 . 

Let us finally return to the Myers-Perry metric (13) which is obtained from the 
5-dimensional Schwarzschild metric 15 and could be relevant for the geometry near 
the horizon of a small r < l braneworld black hole. The cosmological evolution of 
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such black holes formed in the early universe having the metric 


dS\ 





dr 2 + r 2 (dn 2 ) 


(104) 


has been described in details in section 4. It is possible for such black holes to survive 
through the intermediate eras of the universe 20,21 . It is feasible that such black holes 
with masses in the sub-lunar range exist as dark matter in galactic haloes 24,23 . The 
photon sphere for a Myers-Perry braneworlcl black hole in units of 2 M /M 2 is given 

by 


f p h = 



(105) 


and is plotted as a function of black hole mass in Figure 2. The minimum impact 
parameter is 


u m = V2r ph (106) 

The coefficients a and b given by 28 

a = -L (107) 

b = -7T + ^103(4^2) (108) 

can be then be used to construct the strong field deflection angle. Eiroa 28 has 
calculated the positions and magnifications of the relativistic images and compared 
them with those of Schwarzschild black holes. 

Gravitational lensing in the weak field limit by the black hole (104) has been 
worked out 2 '. When the impact parameter u exceeds a few times the horizon radius 
given by Eq.(28), but is still lesser than l (u < l), the application of weak field 



Fig. 2. The photon sphere r p h (in units of the Schwarzschild radius r s ) is plotted versus mass for 
the Schwarzschild black hole (I), and the Myers-Perry braneworld black holes with l = 10 30 /4 (II), 
and l = 10 20 Z 4 (III). 
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lensing could be relevant. The weak field limit of the deflection angle was calculated 
to be 27 


_ 2 MIU 

Ol-weak “TT o' 

M 4 r$ 


(109) 


In order to satisfy the requirement that u < l, and also obtain non-negligible mag¬ 
nification at the image location, the mass of the black hole should be such that 27 


M l 

M 4 l 4 


( 110 ) 


Using the maximum allowed value for l by present experiments 8 , one then obtains 
that M < 10 ~ 8 Mq. So weak field gravitational lensing for such a black hole could be 
applicable only for masses in the sub-lunar range. As discussed above the braneworld 
scenario is conducive to the existence of primordial black holes in the sublunar mass 
range 20,21 . If such black holes exist in our galactic halo, then the magnification 
of their weak field images turns out to be diminished compared to the standard 
Schwarzschild black holes of similar mass 27 . 


Table 2. Estimates for the strong field lensing angle coefficients and ob¬ 
servables for a black hole in the galactic halo with mass M = 10 25 gm 
and D,i = 10 22 cm, for the standard Schwarzschild and the Myers-Perry 
geometry with two different values of l. ( r m = '2.5logTZ.) 


Observables 

Myers-Perry metric 
with extra dimention l 

Schwarzschild 

metric 

l = 1(P U Z4 

l = 10^ u / 4 

9oo(p arc sec) 

0.039 x ltr 12 

0.0123 x 10” 16 

0.03 x 10~ 12 

S(p arc sec) 

0.2042 x 10- 17 

0.0644 x 10- 21 

.0375 x 10“ 15 

r m (magnitudes) 

9.64 

9.64 

6.82 

Um/r s 

3.37 

1.065 x 10 -4 

2.6 

a 

0.707 

0.707 

1 

b 

-0.689 

-0.689 

-0.4002 


Before concluding this section we present a comparitive study of lensing in the 
Myers-Perry geometry (104) with that in the standard Schwarzschild geometry. The 
respective photon spheres are plotted as a function of mass in Figure 2 choosing two 
values of l for the braneworld case. The mass range chosen is the one which could be 
relevant for Myers-Perry braneworld black holes (104) that have been conjectured 
to exist in the form of binaries in the galactic halo 23,24 . Choosing a specific value 
of mass in the above range, the strong field lensing quantities have been evaluated 
separately for the above two metrics. The strong lensing angle coefficients a and 
b. and the observables #oo> 77 and S are displayed in Table 2. It is interesting to 
note that a larger value for the scale of the extra dimension l takes the values of 
the observables #oo and S for the braneworld case closer to those of Schwarzschild 
lensing. This happens because the size of the braneworld black hole which is much 
smaller compared to the Schwarzschild black hole of same mass increases with I for 
fixed mass (28). Of course, the values of these lensing quantities are far beyond the 
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possiblity of verification by present observational capabilities. But the comparison of 
these numbers provides an in principle method of discriminating between different 
gravity models. 


6. Summary and Conclusions 

Braneworld black holes are the potential testing arenas of a rich theoretical structure 
associated with modified braneworld gravity and extra dimensions. If the fundamen¬ 
tal scale of gravity is much lower than the 4-dimensional Planck scale, then it is 
feasible for higher dimensional black holes to be formed in low energy processes. 
A lot of the present interest in the prospect of obtaining observable signatures of 
braneworld gravity is via the properties of the evaporation products of mini black 
holes produced either in particle collisions inside accelerators 9,10 , or in high energy 
cosmic ray showers 11 . Another issue of interest is regarding primordial black holes 
which may be formed through the collapse of overdense regions in the braneworld 
high energy phase of the early universe. Such black holes could have diverse cos¬ 
mological implications. Further, even larger black holes could form by gravitational 
collapse of matter on the brane. Braneworld black holes have entirely different met¬ 
rics compared to 4-dimensional black holes, and carry the signature of the extra 
dimension in their geometries. The focus of this review has been to discuss the cos¬ 
mological consequences of primordial braneworld black holes, and also to analyse 
the features associated with gravitational lensing in several braneworld black hole 
metrics. 

The physics of gravitational collapse on the brane is not yet understood in 
totality 13 . A serious conceptual difficulty arises due to the fact that the gravita¬ 
tional field equations are rendered incomplete on the brane by the effect of the bulk 
gravitational modes. No complete solution to the bulk gravitational field equations 
representing a spherically symmetric vacuum black hole metric on the brane has 
been found till date. Projections of the 5-dimensional Weyl tensor have been used 
in the 4-dimensional brane field equations in various configurations to obtain pos¬ 
sible braneworld black hole metrics. In section 2 we have discussed some possible 
candidate geometries for braneworld black holes. These metrics are in general rather 
different from the standard Schwarzschild metric, and incorporate modifications of 
the standard 1/r gravitational potential in braneworld gravity 6 . The braneworld 
black hole solutions contain interesting features related to the existence of horizons, 
and the modified black hole entropy and temperature. The Reissner-Nordstrom type 
solution 14 contains a negative tidal charge which can be viewed as the reflection of 
the black hole mass by the bulk with negative cosmological constant. Other solu¬ 
tions discussed include the one that incorporates the contribution of the projected 
Weyl tensor in terms of the post-Newtonian parameters 36 . 

The 5-dimensional character of braneworld gravity at small scales motivates the 
consideration of the induced Myers-Perry metric 15 as a braneworld black hole candi¬ 
date for black holes with radius r < l. Physically this corresponds to the fact that a 



February 2, 2008 17:6 WSPC/INSTRUCTION FILE review 


Braneworld Black Holes in Cosmology and Astrophysics 31 

small black hole is unable to distinguish between the bulk dimension and our 3 brane 
dimensions. The modified properties 16 of such black holes within the braneworld 
context have been discussed 47 , and the different subtelties regarding their Hawk¬ 
ing evaporation on the brane and into the bulk have been analysed earler 17 . Special 
features of the interaction of 5-dimensional black holes with the brane have been re¬ 
vealed in the literature 18,25 ’ 48,49,78 . In this review we have highlighted the relevance 
of using the Myers-Perry geometry for the analysis of the processes of Hawking 
evaporation and the accretion of radiation in the early cosmological evolution of 
primordial black holes in the braneworld scenario. 

In section 3 we have provided a short summary of the main features of the 
cosmology in the RS-II model 4 . Essentially the cosmology in the early braneworld 
era (high energy brane regime) is altered by the presence of a term in the Friedmann 
equation on the brane that is quadratic in the energy momentum tensor. This 
modifies the Hubble expansion at times earlier than the time t c = 1/2. Among 
other consequences, such a scenario allows for the possibility of inflation with steep 
potentials that could be excluded in the standard scenario. Various features of 
the cosmology of braneworld models have been analysed in details in some recent 
reviews 1,5 . In the context of the present review, the most significant implication 
of the modified high energy behaviour is that the expansion of the Hubble volume 
makes it feasible for accretion by the primordial black holes of the surrounding 
radiation to take place. This feature is primarily responsible for the mass growth and 
prolonged survival 20,21 of the braneworld black holes. The standard cosmological 
expansion of the universe is recovered at later times, in order for the observationally 
established processes such as nucleosynthesis to work out. 

The description of the braneworld cosmological evolution with primordial black 
holes has been provided in section 4. The modified mass-radius relationship for 
the induced 4-dimensional Myers-Perry black holes leads to slower evaporation and 
longer lifetimes compared to standard Schwarzschild black holes 19 . The detailed 
properties of evaporation of such black holes are described in a recent review 1 '. 
The accretion of radiation in the high energy phase could be effective because the 
growth of black hole mass is smaller than the growth of the mass in the Hubble 
volume in braneworld evolution 20,21 . Thus, a large fraction of primordial black holes 
may survive up to much later eras. The decaying black holes affect several cosmolog¬ 
ical processes at different eras, and hence the observational abundance of different 
species, for example, the background high enery photons 65 , could be used to put 
constraints on the initial mass spectrum of the black holes. These contraints have 
to be evaluated considering the altered cosmological evolution in the braneworld 
scenario, and therefore could be significantly modified compared to similar con¬ 
straints in standard cosmology 22 . The exchange of energy of the black holes with 
the surrounding radiation in the high energy era leads to mass disequilibration 
of neighbouring black holes 23 . As a consequence, binaries could be formed later 
through 3-body gravitational interactions. Such binaries have masses in the sub¬ 
lunar range, and gravitational waves emitted during their coalescence come in the 
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detectable range 24 . 

Braneworlcl black holes exisiting at present times offer another scope of detec¬ 
tion, viz. through the gravitational lensing of light sources by them. In section 5 
we have first briefly reviewed the framework of gravitational lensing using which 
both weak field and strong field lensing can be handled in a unified manner 71 . 
The expressions for the lensing quantities in different braneworld metrics have been 
presented analysing the crucial differences from standard Schwarzschild lensing. 
Myers-Perry black holes existing in the galactic halo in certain mass ranges would 
be difficult to detect via weak field lensing due to reduced magnification compared 
to Schwarzschild black holes 27 . Strong field gravitational lensing through its promi¬ 
nent features such as the production of relativistic images and retrolensing offers the 
ultimate scope of discriminating between different gravity models 28 . Using the black 
hole at the centre of our galaxy as a candidate lense, it is possible to compute the 
theoretical values of several lensing observables in different geometries 29 . We have 
presented the computed lensing observables for a few braneworld metrics. A com¬ 
parison with the lensing observables of the Schwarzschild metric shows that it might 
be possible to distinguish the braneworld black hole with tidal charge 14 from the 
Schwarzschild black hole by the possible measurement of one of these observables 
in the near future. 

The physics of extra dimensions in no longer a field of mere theoretical con¬ 
structs. There is currently a lot of ongoing activity on finding observational signa¬ 
tures of braneworld black holes through different physical mechanisms. In addition 
to the motion of light rays, the properties of massive orbiting particles in these 
geometries are also expected to exhibit interesting features. In particular, it has 
been shown that no stable circular orbits exist in the equatorial plane of the Myers- 
Perry metric 78 . The position of the innermost stable orbit shifts for braneworld 
black holes as is the case with the photon sphere, and this could lead to observable 
modifications of the properties of the accretion disks around black holes. Grav¬ 
ity wave spectroscopy offers another possibility of observing braneworld gravity. 
It has been pointed out using the black string between two branes as a model 
of a braneworld black hole, that the massive bulk gravity modes could lead to 
detectable spectroscopic signatures that are absent in the normal 4-dimensional 
gravitational waves 79 . Furthermore, the possibility of even distingushing between 
different braneworld models such as the ADD model 2 and the RS model 4 via the 
production of their respective black holes in accelerators has been argued 80 . Much 
excitement exists indeed in the prospects of detecting signatures of extra dimen¬ 
sions. The range of cosmological and astrophysical implications of extra dimensions, 
braneworld gravity, and black holes have only started to being investigated. 
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